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We explore the effects of a microscopic nuclear three-body force on the threshold 
baryon density for kaon condensation in chemical equilibrium neutron star mat- 
ter and on the composition of the kaon condensed phase in the framework of the 
Brueckner-Hartree-Fock approach. Our results show that the nuclear three-body 
force affects strongly the high-density behavior of nuclear symmetry energy and 
consequently reduces considerably the critical density for kaon condensation pro- 
vided that the proton strangeness content is not very large. The dependence of the 
threshold density on the symmetry energy becomes weaker as the proton strangeness 
content increases. The kaon condensed phase of neutron star matter turns out to be 
proton-rich instead of neutron-rich. The three-body force has an important influence 
on the composition of the kaon condensed phase. Inclusion of the three-body force 
contribution in the nuclear symmetry energy results in a significant reduction of the 
proton and kaon fractions in the kaon condensed phase which is more proton-rich in 
the case of no three-body force. Our results are compared to other theoretical pre- 
dictions by adopting different models for the nuclear symmetry energy. The possible 
implications of our results for the neutron star structure are also briefly discussed. 
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I. INTRODUCTION 

As well known that around p ~ p where po — 0.16 fm -3 is the empirical saturation 
density of nuclear matter, neutron star matter in chemical equilibrium consists mainly of 
nucleons and leptons, whereas at high enough densities it may become more or less exotic [jj. 
Among the possible exotic phases in dense nuclear matter, the kaon condensation is a subject 
of great interest in ouclear Physics , nadronic physlCS and _ star p^fl B Q. 
Since been proposed by Kaplan and Nelson pj, the possibility of kaon condensation and its 
implications for astrophysical phenomena of neutron stars have been extensively discussed 

by many physicists S 0, 0, 0, . It has been suggested in Ref.Q 

that the condensation of kaons may be understood as a chiral rotation away from a V- 
spin scalar ground state and it is related to a partial restoration of the chiral symmetry 
explicitly broken in the vacuum. The presence of kaon condensation may have important 
consequences for determining the structures and evolutions of neutron stars. For example, 
the kaon condensation may soften substantially the equation of state (EOS) of neutron star 
matter and consequently lower the predicted maximum mass of neutron starsj2|. The phase 
transition from the normal matter to the kaon condensed matter is also expected to affect 
the transport properties and the glitch behavior of pulsars jjj. 

It is shown [5] that the kaon-nucleon sigma term T>xn provides a strongly attractive 
interaction between kaons and nucleons. This attraction reduces the K~ energy u>k in 
nuclear medium and it becomes strong enough at a high enough baryon density above which 
the kaon condensed phase is energetically favorable. It is shown that the kaon-nucleon 
sigma term plays an essential role in determining the formation of such a phase. Besides the 
kaon-nucleon interaction, the high-density behavior of the nuclear symmetry energy is also 
important for determining the threshold density for kaon condensation and the composition 
of the kaon condensed phase. Associated to the high-density symmetry energy there is a 
large uncertainty due to the lack of experimental constraints |18| . The nuclear symmetry 
energy plays its role somewhat in a different way from the kaon-nucleon interaction. It 
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equilibrium neutron star 



21 1 . The electron chemical 



determines to a large extent the electron chemical potentia. 
matter below the critical density for kaon-condensation jigL 
potential p e serves as the energy threshold for kaon condensation since the kaon condensation 
becomes energetically favorable as soon as \i e > ujk- 
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Since realistic nucleon-nucleon (NN) interactions have not been fully determined from the 
chiral theory up to now|22j, various theoretical models of the symmetry energy have been 
adopted for studying the properties of kaon condensation. In the work of Ref . [9] , the simple 
parametrizations proposed by Prakash et al. [3| has been chosen for the nuclear symmetry 
energy. In the investigations of Ref.Q, the symmetry energy derived from the variational 
approach[24| has been applied. It turns out that besides the dependence on T, KN , the 
condensation is also sensitive to the high-density behavior of the symmetry energy. Within 
the framework of a non-relativistic microscopic model based on realistic NN interactions, 
nuclear three-bodv forces are critical for reproducing the empirical saturation properties 
of nuclear matter^ Q Q U Q In onr prevrons worlfi Q, the EOS of nuclear 
matter has been explored in the Brueckner-Hartree-Fock (BHF) approach by adopting a 
microscopic three-body force (TBF) from the meson-exchange current method 



The 

TBF turns out to affect strongly the high-density behavior of the symmetry energy, i.e., it 
makes the density dependence of the symmetry energy much stiffer as compared to the result 
without the TBF contribution. The aim of the present work is to investigate the critical 
density for kaon condensation in neutron stars and the composition of the kaon condensed 
phase in neutron star matter by using the BHF approach for the nuclear symmetry energy. 
In the calculations, we adopt the chiral Lagrangian[5] to extract the kaon-nucleon part of the 
interactions as in Refs-P. lla]. Special attention has been paid on the effects of the nuclear 
three-body force. 

This paper is organized as follows. In Sect. 2 we review briefly the theoretical models 
adopted in our calculations, including the self-consistent BHF approach, the microscopic 
TBF, and the chiral model for kaon-nucleon interaction. Our numerical results are presented 
and discussed in Sect. 3. In Sect. 4 a summary of the present work is given. 



II. THEORETICAL MODELS 

In the kaon condensed phase of neutron star matter, the energy density consists of three 
parts of contributions, i.e. 

£ = £nn + £lep + £kn (1) 

where Snn is the nuclear part, ei ep denotes the contribution of leptons, and €kn is the 
contribution from the kaon-nucleon interaction. 
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Our calculations of the nuclear part of the EOS are based on the BHF approach for 
asymmetric nuclear matter [2^ . The starting point of the BHF approach is the interaction 
G matrix which satisfies the following Bethe-Goldstone (BG) equation [2^ 

nt a \ , \ k ^)Q{k 1 ,k 2 ){kxk 2 \ n( n \ / n 
G{p,P;u) =v NN + v NN ~t\ r G(p,p;u), (2) 

where uo is the starting energy and Q{k%, k 2 ) = [1 — n(&i)][l — n{k,2)\ is the Pauli operator 
which prevents the two intermediate nucleons from being scattered into the states below the 
Fermi sea. The isospin asymmetry parameter is defined as (3 — (p n — p p )/pb, Pn, P P , and 
Pb being the neutron, proton and total baryon number densities, respectively. The single 
particle (s.p.) energy is given by e(k) = h 2 k 2 /(2m) + U(k), where the s.p. potential is 
calculated from the real part of the on-shell G matrix, i.e. 

U(k) = n{k')Re(kk'\G{e{k) + e{k'))\kk') A . (3) 
k' 

In the present calculations, we adopt the continuous choice for U (k) since it has been proved 
to provide a much faster convergency of the hole-line expansion for the energy per nucleon in 
nuclear matter up to high densities than the gap choice (2^. In addition, in the continuous 
choice, the s.p. potential describes physically the nuclear mean field felt by a nucleon in 
nuclear medium. The realistic nucleon- nucleon (NN) interaction vnn is the Argonne Vis 
(AVis) two-body forcej3] supplemented with the contribution of the TBF. 

The microscopic TBF adopted in the present calculations is constructed from the meson- 
exchange current approach [2i|. In this TBF model, four important mesons ir, p, a and lo 
are considered 29 1 . The meson masses in the TBF have been fixed at their physical values 

n 

except for the virtual a- meson mass which has been fixed at 540MeV according to Ref . [28J . 



This value has been checked to satisfactorily reproduce the AV\& interaction from the one- 
boson-exchange potential (OBEP) model 30]. The other parameters of the TBF, i.e., the 
coupling constants and the form factors, have been determined from the OBEP model to 
meet the self-consistent requirement with the adopted AVis two-body force. The values of 
the parameters are given in Ref. [3^ . For a more detailed description of the model we refer 
to Refs. j3, 13- We want to stress that the most important component of this TBF is 
essentially that introduced automatically in the Walecka relativistic mean field theory 
over that in the non relativistic theories. 
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In order to include the TBF contribution into the BG equation, we follow the standard 
scheme 13 to reduce the TBF to an effective two-body interaction which is expressed 

in r space as 

(nr 2 \V.^\r^) = ilr£ J dr^dr^for/Xl - t^ >73 (/ 13 ))(1 - Vr 2 ,r 3 (r' 23 )) 

(4) 

x W 3 (f 1 ' r 2 'r 3 '\fif2f 3 )(p n (T 3 r 3 )(l - r/ ni7S (ri 3 ))(l - i] T2tT3 (r 23 )) 
where the trace is taken with respect to the spin and isospin of the third nucleon and 
the indices ti,t 2 and r 3 denote the isospin ^-components.. The function r) TiTr (r) is the 
defect function. Since the defect function is directly determined by the solution of the 



BG equation 



32|, it must be calculated self-consistently with the G matrix and the s.p. 



potential U(k) [20j at each density and isospin asymmetry. 

By solving self-consistently the coupled Eqs.(J2J), ©, and (@J we can obtain the reaction 
G matrix. From the G matrix, we calculate the nuclear contribution enn to the total energy 
density e for any given baryon density and isospin asymmetry. The Enn of asymmetric 
nuclear matter can be separated into two part, i.e., the isoscalar part and the isovector part. 
The stiffness of the isoscalar part plays an important role in predicting the structure of a 
neutron star while the isovector part is crucial for the chemical properties of neutron star 
matter. Microscopic investigations 20, 30] show that for a fixed baryon density and isospin 
asymmetry, the isovector part of e^n is completely determined by the symmetry energy 
S(u) and can be expressed as p uS(u)(3 2 , where u = Pb/Po is the dimensionless baryon 
density and po — 0.16 fm~ 3 the empirical saturation density of nuclear matter. 

We evaluate the kaon-nucleon sector in the energy density of the kaon-condensed phase 
of matter in terms of the effective chiral Lagrange density which was suggested by Kaplan 
and Nelson 5| and extensively investigated afterwards:?]. In the Kaplan-Nelson model, the 
577(3) x 577(3) chiral Lagrange density is expressed as 
f 2 

C x = —TrdfJJd^U^ + TrB(i^D^ — m B )B 

+ F TrBYldA' B\ + D TrBj^A^ B} 
+cTrM(U + U ] ) + ai TrB{iMi + $Mg)B 

+a 2 TrBB{£M£ + £}M$) + a 3 TrBB Tr(£M£ + (*M&). (5) 

where we adopt the same notations and symbols as in Ref . p] . The above Lagrange includes 
an octet of pseudoscalar masons M. and an octet of baryons B. The first four terms conserve 
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the chiral symmetry. / = 93MeV is the pion decay constant. The constants D and F are 
given by D = 0.81 and F = 0.44. The coefficients aim s ,a2m s , a^,m s , and c determine the 
strength of the chiral symmetry breaking in the last four terms. We adopt a\m s = —67 
MeV and a<ivn s = 134 MeV obtained from the baryon mass splittings 6]. The constant c is 
related to the bare kaon mass by the Gell-Mann-Oakes-Renner relation m 2 K = 2cm s / f 2 . The 
parameter a^m s had remained largely uncertain for many years due to our poor knowledge 
of the strangeness content of the proton and the kaon-nucleon sigma term S m \l . \o . . 
Fortunately, the large ambiguity in this parameter has been settled recently by Dong, Lagae 



and Liu |l6| with small error based on the lattice calculations. In the present calculations 
we adopt the value of a^m s in the range —310 MeV < a^nis < —134 MeV as done in 
Refs.flQ in order to investigate the sensitivity of the TBF effect to the varying of the 
proton strangeness content. Our adopted central value a^rn^ = —222 MeV is very close to the 
value extracted from the lattice guage calculations in Ref. [36] which is —231 MeV with error 
of less than 4%. a^m s measures the proton strangeness content and plays an essential role in 
determining the threshold density for kaon condensation [9] since it provides the dominant 
attraction in the kaon-nucleon interaction via the sigma term Y*kn = — (ai + 2a2 + 4a3)m s /2. 
The above specified range of a^m^ from —310 MeV to —134 MeV corresponds to a reasonable 
range 0.2 > (ss) p / (uu+dd+ss)p> of the proton strangeness content [35j|. Following exactly 
the standard prescript in Ref. |9|, we can get the kaon-nucleon energy density of the kaon 
condensed matter, i.e. 

Skn = sin2 + 2m Kf 2 sin2 f + PB^iXp + 2a 2 + 4a 3 )m s sin 2 ^ (6) 

where ps and x p = p p / Pb denote the baryon number density and the proton fraction, 
respectively. 9 is the amplitude of the condensation and comes from applying the Baym 
th eo rem Q. The lepton part of the energy density £fcp ean be readily obtained in a standard 
way from the contributions of the filled Fermi seas of leptonsp, 16]. 

In the neutron star matter with kaons the chemical equilibrium can be reached through 
the following reactions 

n^p + l + ui, n^p + K~, l^K~ + v u (7) 



where I denotes leptons, i.e., / = e, p. One can determine the ground state by minimizing the 
total energy density with respect to the condensate amplitude 6 keeping all densities fixed. 
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1 equilil 



This minimization together with the chemical e quili brium and charge neutrality conditions 
leads to the following three coupled equations^" 



and 



cos 9 = (^m 2 K f 2 + ^up (2a 1 x+2a 2 +4:a 3 )m s -^fxup (l+x)) , (8) 



9 9 

fx = p e = [Lk = Pn — A%> = 4(1 — 2x)S(u) sec 2 - — 2a\m s tan 2 - , (9) 



f 2 psin 2 9 + up (l + x) sin 2 - - xup + — + t)(\(jl\ - m M ) — = 0, (10) 



A* 3 M n , , W ~ K f' 2 
- - xupo + ~oZ2 W\ ~ m »> ^3 

where p represents the chemical potential of the negative charged particles and u = Pb/Po 
is the baryon number density in unit of po- The last two equations are from the chemical 
equilibrium and charge neutrality conditions, respectively. The EOS and the composition of 
the kaon condensed phase in the chemical equilibrium neutron star matter can be obtained by 
solving the coupled equations (JHJ, ©, and (fTUjl. The critical density for kaon condensation 
is determined as the very point above which a real solution for the coupled equations can 
be found. 



III. RESULTS AND DISCUSSIONS 

The density dependence of the nuclear symmetry energy S(u) is depicted in Fig^ where 
the bold solid line and the thin solid line are obtained from the BHF approach adopting 
the A V\s plus the TBF and the pure A Vis two-body force, respectively. In the figure the 
symmetry energy from other models j^. 1^ are also plotted for comparison. The dashed, 
dotted, and dot-dashed curves correspond to the following parameterizations proposed in 
Ref.|23|, 

S(u) = (2§ - l) ^4 0) (t*l - F(uj) + S F(u) (11) 

and 

2m 2 

F 1 (u)=u, F 2 («) = ^ or F 3 (u) = V^ (12) 

where Sq ~ 30 MeV and = (37r 2 po/2) 2 ^ 3 /2m are the symmetry energy and the Fermi 
energy at the nuclear saturation density, respectively. The double-dot-dashed line is the 
result of Ref.j^J by using the variational approach with the UV14+TNI interaction. It 
is seen from the figure that in both cases with and without the TBF contribution, our 
calculated S(u) are monotonically increasing functions of the baryon density. Below and 
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around the empirical saturation density, the TBF effect is quite small, while at high densities 
the TBF contribution leads to a strong enhancement of the increasing rate of S(u). This 
strong enhancement of the symmetry energy due to the TBF could be explained as follows. 
The symmetry energy (the isovector sector of the EOS) describes the energy required to 
increase the isospin asymmetry of the matter. A higher value of asymmetry implies a larger 
neutron excess in the matter. Therefore for a fixed density, at a higher isospin asymmetry 
the neutron Fermi momentum becomes larger and consequently some nucleons may have 
higher momentum than in the symmetric nuclear matter. Since the TBF is a short-range 
interaction and its effect is stronger for nucleons with lar ger momenta, its contribution to 
the EOS increases as the isospin asymmetry increases j^. 1.3 (J|. 

The density dependence of the chemical potential \x of the negative charged particles in 
/^-equilibrium neutron star matter is given in FigJ2]for three values of a^m s = —310, —222, 
and —134 MeV denoted by a, b, c, respectively. In the figure the solid and dashed curves 
are the results with and without including the nuclear TBF contribution. Below the critical 
density for kaon condensation, the matter is made up of neutrons, protons, and leptons. The 
chemical potential /x is determined by the density dependence of the symmetry energy and it 
increases gradually with density in the case that the symmetry energy is an monotonically 
increasing function of density. The TBF enhances the increasing rate of fi. Above the 
critical density, \i becomes a decreasing function of density. The decreasing rate of /i depends 
strongly on the choice of the proton strangeness content (i.e., on the value of a^m s ), but 
completely insensitive to the choice of different models for the symmetry energy S(u) in the 
case that S(u) is an monotonically increasing function of density as shown in Ref.fl. As 
a consequence, the TBF contribution has almost no any effect on the chemical potential in 
the condensed phase. 

The values of the critical density u c for different models of the symmetry energy S(u), 
and some typical values of the proton strangeness are presented in TablU The results in 
Ref.0] (denoted with F1,F2, and F3) and Ref.Q (denoted with UVU + TNI) are also 
given. As shown in Tabfll the critical density u c for kaon condensation depends on both 
the proton strangeness content (i.e., the value of a^m s ) and the high-density behavior of 
the symmetry energy. The critical density u c is more sensitive to the proton strangeness 
than to the symmetry energy. The value of u c turns out to be lowered by increasing the 
proton strangeness, since the attractive interaction provided by the kaon-nucleon sigma 
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TABLE I: Critical density u c for kaon condensation in unit of pq. 





a 3 m s [MeV] 


Fl 


F2 


F3 


BHF 


BHF+TBF 


UV14+TNI 


a 


-310 


2.4 


2.3 


2.6 


2.6 


2.4 


2.8 


b 


-222 


3.1 


2.9 


3.4 


3.4 


2.9 


4.1 


c 


-134 


4.2 


3.8 


4.9 


5.0 


3.8 





term becomes smaller as the strangeness increases. This is in good agreement with the 
previous investigations by adopting different models for the symmetry energy 

HQ. 

It is 

seen from the table that the predicted values of u c become less sensitive to different models 
for the symmetry energy as the strangeness increases. The TBF affects the critical density 
u c via its contribution to the symmetry energy and its effect is to reduce the critical de nsity 
since it provides an additional repulsive contribution to the isospin symmetry energy [30]. 
Our predicted critical density is in the range of 2.4p — 3.8po if the TBF is included and 
2.6po — 5.0po hi the case without the TBF contribution. In the case of no strangeness, 
inclusion of the TBF contribution in the symmetry energy reduces u c by more than 20% 
from 5 to 3.8. However, if the strangeness is higher, the TBF effect on u c becomes somewhat 
smaller. For instance, in the case of a 3 m s = — 310MeV, the TBF leads to only a less than 
10% reduction of u c . This can be readily understood since on the one hand, for larger 
values of the strangeness, the kaon condensation sets in at lower densities where the TBF 
contribution to the symmetry energy is relatively small. On the other hand, if the strangeness 
is higher, the in-medium energy of K~ drops down faster as a function of baryon density. 
As a consequence, the role played by the strangeness becomes more predominant over that 
by the symmetry energy for a higher value of the strangeness. 

In Fig|31 is plotted the predicted composition of /^-equilibrium neutron star matter for 
three different values of a 3 m 3 = —310 MeV(the upper panel), —222 MeV(the middle panel), 
and —134 MeV(the lower panel) for both cases with the TBF (curves with symbols) and 
without the TBF (curves without symbols). In the figure the proton fraction x p = p p /pb, the 
kaon fraction xk = Pk/pb, and the lepton fraction xi ep = x e + x^ = (p e + p^) / Pb ar e given by 
the solid curves, dashed curves, and dotted curves, respectively. Below the critical density, 
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the matter is in its normal phase which is a highly neutron-rich and charge-neutral mixture 
of neutrons (n), protons (p), electrons (e), and muons (/i). In the normal n,p,e,fi matter, 
the proton fraction x p increases gradually with density and the TBF contribution makes the 
x p rise faster. However, in the condensed phase, the matter becomes more or less proton-rich 
instead of neutron-rich, since a large fraction of protons is required to balance the negative 
charge of K~ which is so abundant in the kaon condensed phase. At high enough density, 
the matter contains even positrons in chemical equilibrium to ensure the charge neutrality. 
These results are in agreement with the previous investigations P. \\\. From FigEl one can 
see that the proton and kaon fractions in the condensed phase depend sensitively on the 
high-density behavior of the symmetry energy. Inclusion of the TBF contribution makes 
the kaon condensed neutron star matter more symmetric in protons and neutrons. As 
the matter goes from the normal phase to the kaon condensed phase, the redistribution of 
charge depends strongly on the high-density behavior of the symmetry energy. If the TBF 
contribution is included, the symmetry energy rises much faster as a function of density. As 
a consequence, the matter is driven to more symmetric in neutrons and protons in order to 
reduce the additional repulsive isospin energy due to the TBF contribution. 

In Ref. [n| the kaon condensation has been studied by adopting the symmetry energy 



from the variational calculations in Ref. 



24j . In contrast to the present results from the BHF 



3811 and the relativistic 



approach as well the predictions from the Dirac-Brueckner method 

mean field theory j^j], the symmetry energy calculated in Ref. 24| decreases with density 
at high enough densities. As a consequence, the composition in the kaon condensed phase 
in Ref . [r| strongly deviates from our result due to the completely different high-density 
behavior of the symmetry energy. For example, in the most extreme case of the UV14+TNI 
interaction, the obtained neutron star becomes almost a "proton" star just above the critical 
density for kaon condensation as shown in Fig. 4. While in the present calculations, the 
neutron star matter with kaons prefers to be symmetric nuclear matter in agreement with 

L 

the results reported in Ref. The discrepancy between the high-density behavior of the 
symmetry energy obtained from the BHF and the variational approaches is still not clear 0, 



271 ] and deserves further investigations. 

It should be stressed that in the present investigations, we do not consider the mixed phase 
of the normal matter and the condensed matter. It has been pointed out [8] that in the first- 
order phase transition to the kaon condensed phase, the conservation law(s) can only be 
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fulfilled in a global rather than a local sense. This means that in the mixed phase, the local 
charge-neutrality, i.e. Eq. (fTU|) becomes invalid. Since our aim is to investigate the nuclear 
TBF effects on the critical density and the chemical composition of the condensed phase, 
the present simplification is somewhat desirable for our purpose. However, for practical 
applications in neutron stars we should consider the problem of the mixed phase formation 
as done in Ref.jsJ which will be discussed elsewhere. 

Before summary, we discuss briefly the possible implications of our results for the maxi- 
mum mass of neutron stars. In FigJS]we display the predicted EOS of neutron star matter 
for both cases with kaons (solid curves) and without kaons (dashed curves). In the figure the 
results by using the AVis plus the TBF (right panel) are compared to the ones by adopting 
the pure AVis two-body force (left panel). It is seen that in the case of no kaons, inclusion 
of the TBF contribution makes the EOS at high densities become much more stiffer. As a 
consequence, the maximum mass of neutron stars becomes much larger if the TBF contribu- 
tion is included as shown in Ref.j3] where the calculated maximum mass is about 2.3 M 
and 1.6 M respectively for the two cases with and without including the TBF in the cal- 
culations for our central value of a^m s = —222 MeV. The latter value is close to the 1.5 M Q 
found in Ref.Jsj for a compression modulus of 210 MeV. In the case that kaons are allowed, 
the EOS of the matter in the condensed phase is softened considerably due to the attractive 
interaction between nucleons and kaons. Since a softer EOS implies a smaller maximum 
neutron star mass, the kaon condensation may lower the predicted maximum mass in both 
cases with and without the TBF. In our present work, however, we do not consider the role 
of other possible strange particles, such as A, E, H hyperons. The appearance of hyperons 
may provide an additional softening of the neutron star EOS and consequently reduce fur- 
ther the predicted maximum mass 41j . One would have to minimize the free energy at each 
density in order to choose the appropriate configuration. 



IV. SUMMARY AND CONCLUSION 



In summary, we have investigated the nuclear TBF effects on the kaon condensation in 
neutron star matter in the framework of the BHF approach. The TBF affects the critical 
density for the kaon condensation through its repulsive contribution to the symmetry energy. 
In both cases with and without the TBF, the calculated symmetry energy is an monotonically 
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increasing function of baryon density, in agreement with the results from the relativistic 
mean field approach and the DB approach. The TBF repulsion turns out to increases 
rapidly as the baryon density increases. As a consequence the high-density behavior of the 
symmetry energy becomes much stiffer as compared to the one obtained by adopting only 
the two-body nuclear force. Our results show that inclusion of the TBF contribution in the 
symmetry energy reduces the critical density by about 8% ~ 25% depending on the choice 
of the proton strangeness content. The influence of the TBF on the critical density becomes 
smaller if the strangeness is higher. The predicted critical density is in the range from 2.4p 
to 3.8p when the TBF is included. 

The additional repulsive contribution to the isospin energy due to the TBF drives the 
kaon condensed phase of neutron star matter to become more symmetric in neutrons and 
protons as compared to the results without the TBF. In the normal phase of neutron star 
matter the proton fraction is small and the matter is highly neutron-rich. The proton 
fraction gradually increases with baryon density. While in the kaon condensed phase, the 
matter becomes proton-rich in order to balance the negative charge of the kaon field. The 
composition in the kaon condensed phase depends sensitively on the high-density behavior 
of the symmetry energy. We find that the additional repulsion from the TBF lowers the 
proton and kaon fractions in the kaon condensed phase and it is so strong at high densities 
as to make the condensed matter almost symmetric in neutrons and protons. The EOS of 
neuron star matter is found to be softened considerably by the kaon-nucleon interaction in 
the kaon condensed phase. 
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FIG. 1: The density dependence of nuclear symmetry energy from different models: Fl(dashed 
curve), F2(dotted curve), F3 (dot-dashed curve), UV14+TNI(double-dot-dashed line), BHF(thin 
solid line), and BHF+TBF(bold solid line). 

FIG. 2: The chemical potential of the negative charged particles for three different values of 
a^m s = —310, —222, and —134 MeV denoted by a,b and c respectively. The solid curves represent 
the results including the TBF contribution, while the dashed curves are the results without the 
TBF contribution. 

FIG. 3: The composition of chemical equilibrium neutron star matter for the two cases with (curves 
with symbols) and without (curves without symbols) the TBF contribution. The solid curves are 
the results for the proton fraction x p , the dashed curves for kaon fraction xk, and the dotted curves 
for lepton fraction xi. 

FIG. 4: The composition of neutron star matter by using different models for the symmetry energy 
and adopting a^,m s = —222 MeV. The curves with symbols are obtained from the BHF calculations 
including the TBF contribution (BHF+TBF). The bold curves without symbols correspond to the 
results of Ref.Jlfl by adopting the symmetry energy of the variational approach and the UV14+TNI 
interaction 



-Q. 



241 ]. The thin curves are obtained by using the linear density dependent symmetry 
energy 9]. The solid and dashed curves are the results for the proton fraction x p , and the lepton 
fraction xi respectively. 



FIG. 5: The predicted EOS of neutron star matter for both cases with kaons (solid curves) and 
without kaons (dashed curves). The results by using the AV\% plus the TBF (right panel) are 
compared to the ones by adopting the pure A Vis two-body force (left panel). 
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